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Abstract 

The performance of the method of angular moments on the Ar s determination 
from analysis of untagged decays B®(t), B s (t) — ► J/ip{-^ K + K~) is 

examined. The results of Monte Carlo studies with evaluation of measurement 
errors are presented. The method of angular moments gives stable results for the 
estimate of Ar s and is found to be an efficient and flexible tool for the quantitative 
investigation of the B® — > J/ip^ decay. The statistical error of the ratio Ar s /r s 
for values of this ratio in the interval [0.03, 0.3] was found to be independent on 
this value, being 0.015 for 10 5 events. 



1 Introduction 

The study of decays B° s (t),lf s (t) -> J/V>H /+r ) K + K~), which is one of the gold 
plated channels for I?-physics studies at the LHC, looks very interesting from the physics 
point of view. It presents several advantages related to the dynamics of these decays, 
characterized by proper-time-dependent angular distributions, which can be described 
in terms of bilinear combinations of transversity amplitudes. Their time evolution in- 
volves, besides the values of two transversity amplitudes at the proper time t = and 
their relative strong phases, the following fundamental parameters: the difference and 
average value of decay rates of heavy and light mass eigenstates of B® meson, Ar s and 
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T s , respectively, their mass difference AM S , and the CP-violating weak phase <f>(*\ The 
angular analysis of the decays B®(t),B~° s (t) — > J/ip(-^ / + Z~)0(— > K + K~) provides com- 
plete determination of the transversity amplitudes and, in principle, gives the access to 
all these parameters. 

In the present paper we examine the performance of the angular-moments method 
[T] applied to the angular analysis of untagged decays B®(t), B° s (t) — > J/ip{— > <p{— > 
K + K~) for the determination of Ar s . After giving the physics motivation in Section 2, we 
describe in the next section the method of angular moments based on weighting functions 
introduced in Ref . P . For the case of Ar s determination this method is properly modified 
in Section 4. The SIMUB-package [2] for physics simulation of 5-meson production and 
decays has been used for Monte Carlo studies with two sets of weighting functions. In 
Section 5 we present the results of these studies and concentrate on the evaluation of 
measurement errors and their dependence on statistics. 



2 Phenomenological description of the decays 

B°(t),B°M -» l + l-)4>{^ K + K~ 



The angular distributions for decays B° s {t),B s {i) — > J/ip{-^> l + l~)cj){-^ K + K~) are gov- 
erned by spin-angular correlations (see jSJ-jE]) and involve three physically determined 
angles. In case of the so-called helicity frame jSj, which is used in the present paper, 
these angles are defined as follows (see Fig. [TJ): 

• The z-axis is defined to be the direction of ^-particle in the rest frame of the B®. 
The x-axis is defined as any arbitrary fixed direction in the plane normal to the 
z-axis. The y-axis is then fixed uniquely via y = z x x (right-handed coordinate 
system) . 

• The angles (0;+, Xi+) specify the direction of the l + in the J/ip rest frame while 
{&K+-, Xk+) give the direction of K + in the <fi rest frame. Since the orientation of 
the x-axis is a matter of convention, only the difference \ = xi+ ~ Xk+ of the two 
azimuthal angles is physically meaningful. 

In the most general form the angular distribution for the decay B®(t) — > J/if)(— > 
l + l~) > K + K~) in case of a tagged B° s sample can be expressed as 

d 4 N ta9 (B°) 9 6 

■ - ^Q^.e^x). (l) 



cfcosO/+ dcosO^ + dx dt 327r ~[ 

Here Oi (i = 1, 6) are time-dependent bilinear combinations of the transversity ampli- 
tudes Ao(t), A\\(t) and A±(t) for the weak transition B®(t) — ► J/ip4> [7 j (we treat these 
combinations as observables) : 

d = |A (t)| 2 , G 2 = |A||(t)| 2 , 3 = \A ± (t)\\ 

A = Tm(Af\(t)A ± (tj), O 5 = Re(A* (t)A ll (t)), O e = lm(A*(t)A ± (t)) , (2) 
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Figure 1: Definition of physical angles for description of decays B®(t), B s (t) — > J/ip( 
(f>(— > K + K~) in the helicity frame 

and the ^ are functions of the angles 0/+, 0^+, x on ly 13 : 



9i = 


2cos 2 0x+sin 2 0;+ , 


92 = 


sin 2 0^-+(l — sin 2 0£+cos 2 x) 


93 = 


sin 2 0^-+(l — sin 2 0;+sin 2 x) 


94 = 


— sin 2 0^:+sin 2 0j+sin2x , 


95 = 


^^sin20 i +sin20^+cosx , 
v2 


96 = 


-^sin20z+sin26£;+sinx . 
v2 



(3) 



For the decay B s (t) — > J/i})(— > l + l ) (f>(— > K + K ) in case of a tagged B s sample the 
angular distribution is given by 

#Ntag(B a ) = 9 ^(^.(Q^e^x) (4) 

dcos<d[+dcost!) K +dxdt 32n ~{ 
with the same angular functions gi and 

1 = |A(t)| 2 , 2 = \A\\{t)\\ 3 = \A ± (t)\ 2 , 

4 = Im(4(t)A ± (t)), £> 5 = Re(^(^||W), O e = lm(A* (t)A ± (t)) , (5) 
where A (t), Au(t) and A^(t) are the transversity amplitudes for the transition B° s (t) — > 
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The time dependence of the transversity amplitudes for the transitions B®(t) , B s (t) — > 
J/ip<f) is n °t °f purely exponential form due to the presence of B® — B^ s mixing. This 
mixing arises due to either a mass difference or a decay-width difference between the 
mass eigenstates of the (B® — B s ) system. The time evolution of the state \B°(t)) of an 
initially, i.e. at time t = 0, present B° s meson can be described in general form as follows: 

\B° s (t)) = g + (t)\B° s ) + g-(t)\B° s ) , 9+{t = 0) = 1 , «/_(t = 0) = , 

i.e., the state \B®(t)) at time t is a mixture of the flavor states \B®) and \B° S ) with prob- 
abilities defined by the functions g+(t) and g~(t). In analogous way, the time evolution 
of the state \B s (t)) of an initially present B s meson is described by the relation 

\B° s (t)} = g + (t)\B s }+g-(t)\B s } , g + (t = 0) = 0, g.{t = 0) = 1 . 
Diagonalization of the full Hamiltonian (see jHj for more details) gives 



g + (t)=g-(t)/a 2 , g.(t) = g+(t) . (6) 

Here Hl/h = M L / H — {i/2)Y L / H are eigenvalues of the full Hamiltonian corresponding 
to the masses and total widths of "light" and "heavy" eigenstates | and a is 
a phase factor defining the CP transformation of flavor eigenstates of the neutral B s - 
meson system: CP\B° S ) = a\~B° s ). In the case \a\ ^ 1 the probability for B® to oscillate 
to a B° s is not equal to the probability of a B° s to oscillate to a 5 S °. Such an asymmetry 
in mixing is often referred to as indirect CP violation, which is negligibly small in case 
of the neutral .B-meson system. 

The time evolution of the transversity amplitudes Af(t) (/ = 0, ||, _L) is given by the 
equations 



A f (t) = A f (0) \g + (t) + g_(t)J—$A , A f (t) = A f (0) 
I r] J CP a J 



Vcp a 



■ (7) 



Here t)q P are eigenvalues of CP-operator acting on the transversity components of the 
final state which are eigenstates of CP-operator: 



cp\j/i(>(P) f = T l f CP \j/i><p) f , (/ = 

and ^ is the CP- violating weak phase 0: 



hp 



Vcp ~ 1 ' Vcp — 1 ' Vcp ~ ^ 



4 S) = , ^ = 2[axg(V t * s V tb ) - arg(V^V*)] = -25 7 , 

where 5 is the complex phase in the standard parameterization of the CKM matrix 
elements Vij (i G {u, c, t}, j G {d, s, b}), and 7 is the third angle of the unitarity triangle. 
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The phase <fyg' is very small and vanishes at leading order in the Wolfenstein expan- 
sion. Taking into account higher-order terms in the Wolfenstein parameter A = sin^c = 
0.22 gives a non- vanishing result |l()j : 

= -2X 2 r] = -2A 2 J R 6 sin 7 . 

Here 

Rh 



1 


\v ub \ 


A 


\Vcb\ 



is constrained by present experimental data as Rb = 0.36 ± 0.08 [TT]. Using the experi- 
mental estimate 7 = (59 ±13)° ^2], the following constrain can be obtained for the phase 

0« = -0.03 ± 0.01 . (8) 

According to Eq. J7J) at time t — 0, the transversity amplitudes of B® , B° s — > J/ip(fi 
decays depend on the same observables |A (0)|, |>W0)|, |Aj_(0)| and on the two CP- 
conserving strong phases, Si = ar g [AM0)Aj_(0)] and 5 2 = arp [Aq(0) A j_(0)]. Time- 
reversal invariance of strong interactions forces the form factors parameterizing quark 
currents to be all relatively real and, consequently, naive factorization leads to the fol- 
lowing common properties of the observables: 

lm[A*(0)A ± (0)} = 0, Im[Afj(0)A ± (0)] = 0, Re[A*(0)A|,(0)] = ±\A o (0)A l{ (0)\ . 

Moreover, in the absence of strong final-state interactions, 5\ = it and 62 = 0. 

In the framework of the effective Hamiltonian approach the two body decays, both 
Bg — > J/ip<p an d B% — > J/if)K*, correspond to the transitions b — *> sec with topologies 
of color-suppressed spectator diagrams shown in Fig. |2j Factorizing the hadronic matrix 
elements of the four-quark operators of the effective Hamiltonian into hadronic matrix 
elements of quark currents, the transversity amplitudes |A (0)|, |Am(0)|, |Aj_(0)| of decays 
B®, B° q — * J/i/;V((q, V) G {(s, 0), (d, K*)}) can be expressed in terms of effective Wilson 
coefficient functions, constants of J /if) decay, and form factors of transitions B q — > V 
induced by quark currents [I] . In Table [T] we collect the predictions of Ref . [T] for the 
transversity amplitudes of B° s — > J/ip (B% — *> J/ipK*) calculated with B — > K* form 
factors given by different models ITU ITH] . The B —>■ K* form factors can be related 
to the B — > (p case by using SU(3) flavor symmetry. The most precise polarization 
measurements performed recently in decays B — > J/ipK*: 

\A (0)\ 2 = 0.60 ±0.04, |v4j_(0) | 2 = 0.16 ±0.03 (BaBar , 
|A (0)| 2 = 0.62 ±0.04, \A ± (0)\ 2 = 0.19 ±0.04 (Belle PU) , 

confirm the predictions based on the model [To] . 



5 




q 



Figure 2: Color suppressed diagrams for decays B® — > J/ipV ((g, V) G {(s, (ft), (d, K*)}) 

Table 1: Predictions for — > J/ip<p (hi brackets - for B® — > J/ipK*) observables 
obtained in Ref. pQ for various model estimates of the B — > K* form factors fSJ El 
(the normalization condition |v4 (0)| 2 + |v4||(0)| 2 + |v4j_(0) | 2 = 1 is implied) 



Observable 


BSW [13 


Soares [Hj 


Cheng [15 




A)(0)| 2 


0.55 (0.57) 


0.41 (0.42) 


0.54 (0.56) 




^(0)| 2 


0.09 (0.09) 


0.32 (0.33) 


0.16 (0.16) 



3 Angular-moments method 

The angular distributions for decays B° s (t),B s (t) -> J/^(-> l + l~) K + K~) m case 
of tagged 5° and 5^(t) samples (see Eqs. (JTJ) and (jU), respectively) as well as in case of 
the untagged sample can be expressed in the most general form in terms of observables 
k(t): 

mi + ,0 K+ ,x,t) = ^ EM*)#(©z+: (9) 



32tt 



i=i 



The explicit time dependence of observables is given by the following relations: 



h(t) = 


|^o(0)| 2 G L (t), 


hit) = 


l^||(0)| 2 G L (t), 


hit) = 


|^±(0)| 2 G H (t), 


hit) = 


1^1(0)11^(0)1^), 


hit) = 


\A o {0)\\M°)\ G L(t) cos(«J 2 - <5i) , 


hit) = 


\A (0)\\A ± (0)\Z 2 (t), 


where we have used the general compact notations: 


G L /nit) = - 


"(1 ± cos0( s) )e- r ^ + (1 =f cos0( s) )e 


Z 1>a (t) = i(e 


-r H t _ e - r ^*)cos5i i2 sin^ 



-rvt 



(10) 
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- for observables 6, = (Oi + Cj)/2 in case of the untagged sample with equal initial 
numbers of B® and ~B° S , while 

= G L/H (t) ± e~ r ^ sin(AMt) sin</>« , 
z[ B 2 ° y(S ° s \t) = Z lt2 (t) ± e" rs * [sin5 li2 cos(AMt) - cos5i, 2 sin(AMt)sin^ ) 

- for observables b\ B "^ = Oi and b\ Bs ^ = Oi in case of tagged and B° s (t) samples, 
respectively, with T s = (r L + T^)/2. It is easy to see that both in the tagged and 
untagged case we have 

G ^WI^) =0 = e- r ^. 

According to Ref. [I], the observables bi(t) can be extracted from distribution func- 
tion © by means of weighting functions Wi(Qi+, ©k+, x) f° r eac h i such that 

9 r 

— j dcosQi+dcosQ K +dx Wi(Q l+ ,Q K +, x) #i (©/+,©!<:+, x) = kj , (11) 
projecting out the desired observable alone: 

bi(t) = J dcosQi+dcosQ K +dx /(©;+, @^+, X ! *) ®k+, x) ■ (12) 

The angular-distribution function Q obeys the condition 

L{t) = J dcosG l+ dcosQ K+ dx f(®i+,® K +, X ]t) = M*) + &2(t) + &s(t) • (13) 

For decays i? — > J/ip{-^ Z + /~)0(— > K + K~), the explicit expressions of weighting 
functions, given in Table 5 of Ref. pQ for physically meaningful angles in the transversity 
frame, get the following form (Set A) after transformation into the helicity frame: 

W { A ^ = 2-5 cos 2 B l+ , 

= 2 — 5 sin 2 0;+ cos 2 x > 

(A) 2 r\ 2 

W3 = 2 — 5 sin <di+ sin x > 
= — - sin 2 ®k+ sin 2x , 

= — -1= sin 2Q K + sin 20^+ cos x , 
4y 2 

Wq A ^ = — -1= sin 20^+ sin 20;+ sin x ■ (14) 
4v 2 

The expressions of Eq. (114)1 are not unique and there are many legitimate choices of 
weighting functions. A particular set can be derived by linear combination of angular 
functions g>j (see [T] for more discussions): 

6 

Wi(Qi+, Q K +, x) = J2 X ij9j(<S>i+, <S>k+,x) , (15) 
5=1 



where the 36 unknown coefficients \j are solutions of 36 equations 
9 



J2 \j / dcosQ l+ dcosQ K + d\ g j (Qi+,Q K +, x) gk(<S>i+,<S>K+,x) = 8 ik ■ (16) 
Sin j=1 J 

The weighting functions (set B) corresponding to the linear combination of the angular 
functions (J3J) are given by 

w[ B) = ^[28cos 2 e x+ sin 2 e,+ - 3 sin 2 6^+(l + cos 2 6,+ )] , 

w {B) = --[4cos 2 e x+ sin 2 e /+ - 29 sin 2 6^+ (1 - sin 2 9,+ cos 2 x) 
8 

+ 21 sin 2 e K +(l- sin 2 6,+ sin 2 x)) , 
w ( 3 B) = - - [4 cos 2 Q K + sin 2 9 /+ + 21 sin 2 Q K +(1 - sin 2 G l+ cos 2 x) 



- 29 sin 2 e K +(l - sin 2 6;+ sin 2 X )} , 
w{ B ^ = sin 2 0^+ sin 2 9;+ sin 2% , 



(B) _ (A) 
5 — w b i 



<*> = < ] - (17) 

For a limited number of experimental events N in the time bin around the fixed value 
of the proper time t, distributed according to the angular function (PJ, it is convenient 
to introduce the normalized observables 

bi(t) = bi(t)/L(t) (18) 

with normalization factor L(t) given by Eq. (fT3*|) . Then, as it follows from the Eq. (fT2|) . 
the observables t>i(t) (fTHj) are measured experimentally by 

1 N 

bt P) = ^E«4 (19) 

with summation over events in a time bin around t. Here w\ = Wi(Oj + , Q J K +, X"0> where 
0;+, Q 3 K + and are angles measured in the j-th event. The statistical measurement 
error of the observable (fTTIj) can be estimated as 



6l (exp) = 1 

N' 



N 



with summation over all events in the same time bin. 

4 Time-integrated observables 

For data analysis it is rather convenient to use the time-integrated observables defined as 
1 f T o r 

h(T ) = -~— / dt / dcosQi+ dcosQ K + dxWi(Qi+,e K +,x) f(®i+, @^+, X! *) (20) 
L[T) Jo J 
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with argument T < T, where T is the maximal value of the 5-meson proper time 
measured for the sample of events being used, and L(T) is a new normalization factor, 
which has the form: 



L(T) 



dt / dcosQi+ dcos®K+ d-X f( cos< 3>i+, cos9^+, x, t) 



L(t) 



\A (0)\ 2 + |A|,(0)| 2 )G L (T) + \A ± (0)\ 2 G H (T) 



(21) 



where, in the compact notations used in Eq. pOjl . 



Gl/h(T) 



T 



dtG L/H {t). 



The following normalization condition is valid for the observables (|20|) : bi + b 2 + b% = 1. 
For a limited number of experimental events N(T), measured in the proper time region 
t F [O.T]. Eq. piljl reduces to 



bt P) (T ) 



N(T Q ) 



N(T) § 
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:22i 



with summation over all events N(Tq) in the time interval t G [0,Tq] for To < T. In case 
of the untagged sample we have 



G L/H {T) 



'licos^ 



+ (1 =F cose 



H 



(23) 



and 



Z(T) 



1 

cos^i^ sin 
1 



dtZ l>2 (T) 



i)/r 2 



For the untagged sample the explicit form of time-integrated normalized observables 
(|2~U]l in terms of the functions Gl/h(T) and Z(T) is given by 



6i(?o) 


= |A (0)| 2 G L (T )/T(T), 




& 2 (T ) 


= |An(0)| 2 G L (T )/T(T), 




& 3 (T ) 


= |A ± (0)| 2 Gh(T )/T(T), 




& 4 (T ) 


= |A||(0)| |A ± (0)| Z(T ) cos^i sin^ 


^/L(T) , 


& 5 (T ) 


= |Ao(0)||A||(0)|G L (T )cos(5 2 - 


*0/£(r) , 


& 6 (T ) 


= |A (0)| |A ± (0)|Z(T )cos5 2 sin(; 


>«/Z(T) . 



(24) 

In the Standard Model (SM) sin0^ ~ and the observables 64,5 (T ) are vanishing. In 
case of a new physics signal the values of sin0^ and 64,5 (T ) can be sizable, however. 
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The following relations are valid for the observables 
h(T ) = cosd 1 sin0«Z(T ( 





b 3 (To) 


\ Gl(T ) 


G H {T ) 








h(T ) 


\ G L (T ) 


G H (T ) 



6 6 (To) = cos5 2 sin0WZ(T o 

If we introduce the function 

j(T) = G H (T)/G L (T), (25) 

then, the values of initial transversity amplitudes at t — and the strong-phase difference 
(82 — 81) are determined from the observables hiT) = 6j(T = T ) by 

u ,m,2 bi(T) 



6 1 (T)+6 2 (T)+6 3 (T)/7(T) 



6 1 (T)+6 2 (T)+6 3 (T)/7(T)' 



M MM2 &3(T)/t(T) 



6 1 (T)+6 2 (T)+6 3 (T)/7(T) 



sm<^ > cos(5 li2 = = = ? — 



cos(5 2 - 50 = 5V / , (26) 
s/h(T)b 2 (T) 

where we consider the initial amplitudes normalized as |Ao(0)| 2 + |A||(0)| 2 + |Aj_(0)| 2 = 1. 
We have also: 

I 1 r r\ . ln.(T\ n...CT\ 

(27) 

For extraction of the .B^- width difference Ar s = — from experimental data it 
is convenient to use a special set of the time-integrated normalized observables: 

bi(T ) = -~— / dt dcosQ l+ dcos<d K +dxWi(Oi+,0 K +,x) e r * /(©/+, <3>k+, X, *) , 
L(T) Jo J 

(28) 

where T' is some arbitrary initial approximation of the i?°-meson total decay width. 
These observables can be extracted from the experimental events N(T), measured in the 
proper time region t G [0, T], by using the formula 

b t P \To) = 1 ^ f - ) E Wi, (29) 

where W( = e r ' p wj, and summation is performed over all events iV(T ) in the time 
interval P G [0, T }. 
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For the untagged sample, the explicit expressions for the time-integrated observables 
can be easily obtained by replacing b iy Gl/h and Z in the expressions of Eq. (j23| by 
h, Gl/h and Z, respectively, (with the same normalization factor (J2TJ) after introducing 
the following notations 



G L/H (T) = £ dte r,t G L/H (t) 

e Ar L T/2 _ l Q -Ar H T/2 _ 1 



= (i±oos0W)— ^ (l=Fcos^W) — , (30) 

Z(T) = - - f T dte r,t Z h2 (T) 
cosoi i2 sm0c J0 

1 _ e Ar L T/2 1 _ e -AT H T/2 

where ATl/h are auxiliary parameters given by 

at i = 2(r' - r L ) , Ar„ = -2(r' - r H ) . (31) 

Eq. (J2fij) is also valid after such a replacement. 



5 Monte Carlo studies 

For Monte Carlo studies of the estimation of physical parameters by applying the angular- 
moments method, untagged samples of events of B®(t) — > J/ip4> decays have been gen- 
erated by using the package SIMUB j2] with various sets of the input values of initial 
amplitudes |Aq(0)| and |Aj_(0)| and AT S . Other parameters are fixed as follows: 

6 1 = n, S 2 = , T s = \/t s = 2.278 [mm/c]" 1 , <p[ s) = 0.04 . 

The value of T s used corresponds to the lifetime t s = 1.464 ps [12] while the CP-violating 
weak phase </>^ was fixed as the upper limit of the constrain (jHj). The value of AT S is 
expected to be negative in the SM. The combined experimental result for |Ar s |/r s is not 
precise: |Ar s |/T s < 0.52 at 95% CL [T2]. In the approximation of the equal and B® 
lifetimes, the |Ar s | extraction can be improved ^2] : |^r s |/r s < 0.31 at 95% CL. A set 
of the untagged-event samples has been generated with Ar s /T s e [—0.3, —0.01] to study 
the influence of Ar s value on the estimation of B®(t) — > J/ip<f) decay parameters from 
data analysis. 

The values of the time integrated observables b[ exp \T ), defined by Eq. ()20jl. can be 
extracted from data according to Eq. (|22|) by summation of weighting functions for each 
event. The statistical error of 6j(To) is defined by 



\ (stat) 



N{T)\ 



N(T ) 

E (bt P) ~ , (32) 
i=i 
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while a systematic error due to limited precision of angular measurements can be esti- 
mated as 



(Sbi 



\(sys) 



(33) 



Here 



A 3 



dcosQi 



-A(cos 0JH 



+ 



dw 3 



d cos 9^+ 



A (cos Qr+ ] 



In a similar way, the values of the observables b\ exp \T ), defined by Eq. (}28|) . can 
be extracted from the data according to Eq. (|29jl . The formulae for statistical and 
systematic errors for bt P) (T ) can be obtained by replacement of w\ to W- in Eq. 
and the following redefinition of A 3 in Eq. (}3*3*|) : 



+ 



<9u^ 1 2 



Ai 



9 cos 9/+ 



A (cos 9 



+ 



-A(cos Qr+ 



+ 



dW 3 
dx 



A(x) 



+ 



A(t) 



_<9cos9#-+ 

Eq. can be applied to estimate the systematic errors related both to the mea- 
surement precision of the detector and to the limited resolution of the Monte Carlo 
generator. In the SIMUB generator, for each variable V G {cos0;+, cos0^+, x, t} ran- 
domly generated for decays B®(t) — ► J/ijj{— > l + l~) </)(—>■ K + K~), the number of bins 
in the region [K n j n ,K nax ] was set as N gen = 50 000. The generation precision for the 
variable V is defined as A(V) = (V max — V min ) /N 9en and systematic errors (J3~3*)) are pro- 
portional to (N^)- 1 / 2 . The B -meson proper time was generated within the interval 
t G [0,T = 2mm/c] which includes 99.3% of all S-decays. We have used samples with 
a maximum of 100 000 events of the decay B° s — > J/ipcf) because a statistics of about 
83800 events is expected to be obtained per year at the CMS detector at the LHC low 
luminosity under realistic triggering conditions [TUj . 



Table 121 shows the values of the observables b^ xp ^ (T) = b\ exp> (T = T) extracted 
from the Monte Carlo data by applying the sets A and B of weighting functions, given 
by Eqs. (JUjl and (|17|h respectively. Various theoretical models for estimation of the 
transversity amplitudes |v4 (0)| and |v4j_(0)| (see Tabled)) have been considered to fix these 
parameters in the SIMUB generator. It can be seen from Table[2]that the choice of N 9en = 
50 000 provides negligibly small systematic errors for the observables as compared with the 
statistical ones. Moreover, both errors slightly depend on the values of the observables. 
For observables obtained by using the set-B weighting functions, the statistical errors are 
significantly smaller than in case of the set-A weighting functions. We should also note 
that even with the statistics of 100 000 events, the values of observables b^ p \T) and - 
as consequence of Eq. (J27|) - the combination cos^i^ sin0 cannot be extracted from the 
data if the CP- violating weak phase <fyf> is small according to the SM expectation (JSJ). In 
this case, these parameters can be estimated only by using a statistics which is not less 
than 3 x 10 9 B° s {t) -»■ J/^(p decays. 

Analysis of the same Monte Carlo data leads to similar conclusions concerning the 
behavior of statistical and systematic errors for the observables b\ exp \T) = b\ exp \To = T). 
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Table 2: Comparison of the observables b\ e (T), extracted from Monte Carlo data, with 
their values bf h \T) corresponding to various theoretical models for |A)(0)| and |Aj_(0)|. 
A sample of 100 000 decay events generated with Ar s /r s = —0.15 was used. The first 
errors are statistical while the second errors correspond to the systematic uncertainties 



a) BSW model p 


□|: 


i 


bf h) (T) 


6^ p) (T) (set A) 


b { r p) {T) (set B) 


1 

2 
3 
4 
5 
6 


0.5425 
0.3551 
0.1024 
-0.00055 
-0.4389 
0.00067 


0.5409 ±0.0044 ±0.0003 
0.3619 ±0.0047 ±0.0004 
0.0972 ±0.0049 ±0.0004 
-0.0017 ±0.0037 ±0.0004 
-0.4344 ±0.0050 ±0.0003 
0.0037 ±0.0055 ±0.0003 


0.5432 ±0.0024 ±0.0002 
0.3579 ±0.0036 ±0.0004 
0.0991 ±0.0034 ±0.0004 
-0.0021 ±0.0033 ±0.0003 
-0.4344 ±0.0050 ±0.0003 
0.0037 ±0.0055 ±0.0003 



b) Model by Soares [H]: 



i 


6f ] (T) 


~b [ r P) {T) (set A) 


b [ r P) {T) (set B) 


1 


0.3908 


0.3900 ±0.0046 ±0.0003 


0.3955 ±0.0023 ±0.0002 


2 


0.2574 


0.2617 ±0.0049 ±0.0004 


0.2551 ±0.0037 ±0.0004 


3 


0.3518 


0.3483 ±0.0047 ±0.0004 


0.3509 ±0.0037 ±0.0004 


4 


-0.00086 


-0.0083 ±0.0040 ±0.0004 


-0.0017 ±0.0035 ±0.0003 


5 


-0.3171 


-0.3156 ±0.0052 ±0.0003 


-0.3156 ±0.0052 ±0.0003 


6 


0.0011 


0.0008 ±0.0052 ±0.0003 


0.0008 ±0.0052 ±0.0003 






c) Model by Cheng 




i 


bf l \T) 


Vf xp \T) (set A) 


~b { r P) {T) (set B) 


1 


0.5271 


0.5228 ±0.0045 ±0.0003 


0.5267 ±0.0024 ±0.0001 


2 


0.2928 


0.2980 ±0.0048 ±0.0004 


0.2950 ±0.0036 ±0.0004 


3 


0.1801 


0.1791 ±0.0048 ±0.0004 


0.1778 ±0.0035 ±0.0004 


4 


-0.00066 


-0.0030 ±0.0037 ±0.0003 


-0.0034 ±0.0034 ±0.0003 


5 


-0.3928 


-0.3927 ±0.0051 ±0.0003 


-0.3927 ±0.0051 ±0.0003 


6 


0.00088 


-0.0019 ±0.0054 ±0.0003 


-0.0019 ±0.0054 ±0.0003 



Table 3: Comparison of the values of observables b[ exp \T), with V = T s , extracted from 
the Monte Carlo data, with their values bf h \T) corresponding to the model of Cheng 
[T5] for initial transversity amplitudes 



i 


bf h) (T) 


Vf xp) {T) (set A) 


bf xp) {T) (set B) 


1 


2.2036 


2.176 ±0.044 ±0.003 


2.206 ±0.026 ±0.001 


2 


1.2242 


1.282 ± 0.045 ± 0.004 


1.245 ±0.034 ±0.004 


3 


0.9187 


0.917 ±0.045 ±0.004 


0.930 ±0.034 ±0.004 


4 


-0.0073 


-0.099 ±0.036 ±0.003 


-0.094 ±0.032 ±0.003 


5 


-1.6425 


-1.618 ±0.048 ±0.003 


-1.618 ±0.048 ±0.003 


6 


0.0098 


0.067 ±0.050 ±0.003 


0.067 ±0.050 ±0.003 



To illustrate the performance of our method in this case, only the results obtained for 
transversity amplitudes, corresponding to the Cheng's model [T3] , are shown in Table El 
Fig. El shows the dependence of the observables bi(T) and 

K{T)= l ~^ T k{T) (. = 1,2,3) 
on the value of the ratio Ar s /T s . For Ar s = we have 

&i,2,3COUr s =o = &wCOIat.=o = |A),||,±(0)| 2 . 

The observables bi(T) slightly depend on Ar s . The rather strong dependence of the 
observables h{T) on the decay width difference Ar s , shown in Fig. El can be used for 
extraction of this parameter from the data analysis as it will be discussed below. 
Under the assumption (f>(*> = 0, we have from Eq. (J30j) : 

Gf IB {T) = ±2 e 



L/HK ' ~ AT L/H ■ 

Therefore, the values of the auxiliary parameters AT L / H , defined by Eq. (pTTj). can be 
determined separately by using the ratios of observables b\ e (T)/b\ exp \To), extracted 
from the data analysis, and solving numerically the equations which arise from one of 
the following relations: 

bi(T)/bi{T ) = G L (T)/G L (T Q ) (i = 1, 2, 5) (34) 

- to determine Ar^, and the relation 

h(T)/b 3 (T ) = G H (T)/G H (T ) (35) 

- to determine Ar#. Then, the decay- width parameters r s , Ar s and T^/h can be 
determined via T' and ATl/h as 

Ar L - AT H Ar L + Ar H r'^ ATL (w\ 

T s = r , Ar s = , r L/H = r t ~y~ ■ ( 36 ) 
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Figure 3: Dependence of the observables bi(T) and b[(T) = 6j(T)[l — exp(— TT)]/(TT) 
(i = 1,2,3) on the value of Ar<,/r s . The observables have been calculated for the case 
of Cheng's model for transversity amplitudes. 

So, by using some reasonable approximation for T' as a starting point for the data analysis, 
the experimental value of T s can be essentially improved simultaneously with determi- 
nation of Ar s . The statistical error of T s determination is expected to be twice smaller 
than for Ar s determination. 

The direct numerical calculations have shown that the difference between the values 
of observables bi(T) (i = 1,2,3,5), calculated with (jyg' = and <fi^ = 0.04, does not 
exceed 0.01%. Even in case of statistics of 100 000 events this difference is negligibly 
small as compared with statistical errors for these observables (see Table |3J) • Therefore, 
the assumption <// s ) = is a good approximation for T s , AT S and T^/h determination by 
the method described above. 

Table HI shows the results of determination of the decay-width parameters after ap- 
plying the described procedure to the Monte Carlo data. The sample of 100 000 events 
generated in case of Cheng's model with Ar s /r s = —0.15 has been used. Both sets A and 
B of weighting functions have been applied to extract the observables b[ exp \ The value of 
T', which is treated as some arbitrary initial approximation for the total decay width of 
B^-meson, was fixed as T' = 1.05 T s , i.e. it was shifted by 5% relative to the "true" value 
of T s fixed in the Monte Carlo generator SIMUB. The value of T = 0.1 T was chosen as 
it provides the minimal statistical errors to determine the ratios b\ exp \T)/b[ exp \To). In 
Table |U we present the result for AI\ obtained from the ratio Si (T)/&i (To) only, 
which gives the best precision. Table 0] shows that the set-B weighting functions give 
more precise and stable results than the set-A functions. 

To improve the precision of Ar s determination, the same procedure should be re- 
peated with T' fixed to be equal to the value of T s determined at the first step. Because 
of Ar s = AT L = AT H in case of V = T s , the value of Ar s is defined at the second 
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Table 4: Results of determination of the decay-width parameters (in units (mm/c) -1 ) 
based on extraction of the observables b\ exp ^ from analysis of 100 000 Monte Carlo events. 
The input value of Ar s corresponds to AT S /T S = —0.15 



Parameter 


Input value 


Measurement (set A) 


Measurement (set B) 


Ar L 


-0.1139 


-0.103 ±0.058 ±0.003 


-0.110 ±0.034 ±0.002 


AT H 


-0.5696 


-0.478 ±0.137 ±0.012 


-0.554 ±0.101 ±0.012 


T L 


2.4493 


2.444 ±0.029 ±0.002 


2.447 ±0.017 ±0.001 


T h 


2.1076 


2.154 ±0.068 ±0.006 


2.115 ±0.050 ±0.006 


r s 


2.2784 


2.299 ±0.037 ±0.003 


2.281 ±0.027 ±0.003 


Ar s 


-0.3418 


-0.290 ±0.074 ±0.006 


-0.332 ±0.053 ±0.006 



step to be equal to the value of AIY determined from the ratio }}[ xp \t) /b^ xp \T ) using 
Eq. ([34)1 . Using the values of Ar s from Table 0] as an input value of T' at the second 
step, we have obtained finally the following results (to be compared with the input value 
AT S = -0.3418 set in the SIMUB generator): 

AT ex P = _ .33o ± 0.057 ± 0.004 (set A) , 
AT ex P = _ .338 ± 0.034 ± 0.002 (set B) . 

This way one can reduce not only the statistical error but also essentially improve the 
stability of the AT S result even in case of using the set-A weighting functions. 

Table |5] shows the statistical errors of Ar s /T s determination by the described ap- 
proach applied to different statistics of Monte Carlo events generated with various "true" 
values of Ar s . The lack of numbers in the table corresponds to cases when the approach 
is not able to give a certain result for Ar s . The use of set-B weighting functions gives 
more stable results even in case of too small statistics and values of AT S , for which the 
same approach does not work with set-A functions. The sensitivity of the method is 
measured by the statistical error of Ar s /r s , which only slightly depends on the value of 
this ratio and is proportional to 1/ \^N, where N is the number of events. In particular, 
for a statistics 100 000 events, the statistical error is about 0.015, while for 1000 events - 
about 0.15. 

In principle, the value of AT S can be determined similarly by using the ratios b\ expS> (T) fb\ 
or ty exp ' (T) /b\ exp \T) , extracted from the data analysis with T' = T s , and solving the 
equations arising from the relations 

b~i(T)/bi(T ) = G L (T)/G L (T ) (i = 1,2,5), b s (T)/b 3 (T ) = G H (T)/G H (T ) 

or 

bi(X)/bi(T) = G L (T)/G L (T) (i = 1, 2, 5) , b 3 (T)/b 3 (T) = G H (T)/G H (T) . 



16 



Table 5: Statistical errors of Ar s extraction (in units [mm/c] - ) obtained by applying 
the angular-moments method with set-B (set-A) weighting functions to the Monte Carlo 
data samples with different numbers of events 



Aiyr s 


200 events 


500 events 


10 3 events 


10 4 events 


10 5 events 


-0.03 








0.035 (-) 


0.014 (0.023) 


-0.05 








0.046 (-) 


0.014 (0.022) 


-0.1 






0.11 (-) 


0.046 (0.079) 


0.014 (0.024) 


-0.15 






0.13(0.19) 


0.045 (0.078) 


0.014 (0.024) 


-0.2 




0.23 (-) 


0.12(0.18) 


0.048 (0.072) 


0.015(0.026) 


-0.3 


0.21 (-) 


0.23 (-) 


0.18 (0.20) 


0.050 (0.083) 


0.016(0.028) 



But in both these cases the precision of AT S determination turns out to be worse than in 
the approach based on the ratios bf xp \T) /Sj e (T ) because of the weak AlVdependence 
of the bi(T) observables. 

The initial transversity amplitudes and strong-phase difference can be recalculated 
from the values of observables b[ exp \T) according to Eq. (f2T?j) . The results of such deter- 
mination of the parameters |A/(0)| 2 (/ = 0, ||, _L) and cos (£2 — <$i) are shown in TableElfor 
different statistics. We have used the Monte Carlo sample generated with the theoretical 
values of the amplitudes |A)(0)| and |Aj_(0)| corresponding to Cheng's model [TH] . To 
extract the observables b\ exp \T), the set B of the weighting function has been applied to 
Monte Carlo data. To estimate the statistical errors for parameters |Aj(0)| 2 (/ = 0, ||, _L) 
and cos(#2 — Si), the standard error-propagation method has been applied to the statisti- 
cal errors of the observables uf (T), taking into account the correlation between pairs 
of different observables. The systematic errors of the observables related to the limited 
generator resolution are neglected. The total errors for parameters | ^4/ (0) | 2 (/ = 0, ||, J_) 
should also include the additional uncertainty related to the error of calculation of 7(T) 
caused by the error of Ar s (see definition of j{T) in Eq. (|23|) and Eq. (|2H1))- In Table El 
we also show these errors calculated by assuming Ar s = — 0.15T S (see Table EJ) 

5(AT S ) _ [ 30% for 10 000 events, 

| AT, | ~ \ 9.3% for 100 000 events. ^ ' 

6 Conclusion 

For the decay B° s — > J/ip <ft in the framework of the method of angular moments a non-fit 
scheme for separate estimation of the parameters Ar s , T s and |A/(0)| 2 (/ = 0, ||, _L) has 
been proposed, based on analysis of an untagged sample, and studied by Monte Carlo 
method. A strong dependence of statistical measurement errors on the choice of the 
weighting functions has been demonstrated. The statistical error of the ratio Ar s /T s for 
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Table 6: Determination of initial transversity amplitudes and strong-phase difference by 
using the values of observables b\ exp \T) extracted from Monte Carlo data. The events 
sample has been generated for the case of Cheng's model ^5] for transversity amplitudes 
and with Ar s = 0.15 T s . The first errors are statistical while the second errors are caused 
by uncertainties of Ar s determination 



Parameter 


Input value 


10 000 events 


100 000 events 


|A(0)| 2 


0.54 


0.527 ±0.007 ±0.012 


0.5398 ±0.0023 ±0.0011 


|A|,| 2 


0.30 


0.337 ±0.011 ±0.008 


0.3023 ±0.0036 ±0.0006 




0.16 


0.136 ±0.010 ±0.020 


0.1579 ±0.0032 ±0.0018 


cos(5 2 - #l) 


-1 


-1.021 ±0.044 


-0.9962 ±0.015 



values of this ratio in the interval [0.03, 0.3] was found to be independent on this value 
and about 0.015 for 10 5 events. The method of angular moments gives stable results for 
the estimate of Ar s and is found to be an efficient and flexible tool for the quantitative 
investigation of the B° s — > J '/ 'ijj <j) decay. 

We would like to thank G. Bohm for the careful reading of the manuscript and useful 
discussions. 
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